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We investigate collisions of polar molecules in quasi-2D traps in the presence of an external 
electric field perpendicular to the collision plane. We use the quantum-defect model characterized 
by two dimensionless parameters: y and s. The first of them is related to the probability of 
the reaction at short distances, whereas the latter one defines the wave function phase at short 
distances. For y close to unity we obtain universal collision rates determined by the quantum 
reflection process from the long-range part of the interaction potential that depends only on the 
van der Waals coefficient, dipole-dipole interaction and the trap frequency. For small short-range 
reaction probabilities collision rates are highly nonuniversal and trap induced shape resonances are 
visible. For high dipole moments we observe the damping of reactive collisions, which can stabilize 
the ultracold gas of polar molecules. The calculations were performed with help of multichannel wave 
funcion propagation by imposing short-range boundary condition derived from the quantum-defect 
model. 

PACS numbers: 34.50.Cx,34.50.Lf 


I. INTRODUCTION 

Successful realization of relatively dense ultracold 
gases of molecules in traps w and optical lattices EHH] 
in their rovibrational ground state offers new possibili¬ 
ties for physics and chemistry [MI]. Ultracold molecule 
reaction rates can be relatively large as it was recently 
demonstrated for fermionic "^^K^^Rb [111113]. This highly 
reactive molecule belongs to the class of molecules, that 
have universal elastic and reactive collision rates. They 
can be calculated analytically from the knowledge of the 
long range part of the interaction potential only, pro¬ 
vided the reaction probability at small distances is equal 
to unity [M] [T^. When reaction probability at short 
distances is smaller than one, the collision rates become 
dependent on the short range phase and the effect of res¬ 
onances have to be taken into account [IMl8] . 

Application of an external electric field to a gas of ro¬ 
tationless polar molecules polarizes them, and as a con¬ 
sequence increases the rates of reactive collisions lead¬ 
ing to significant reduction of dipolar gas lifetime [13 
[13 ITB] . The essential condition for creation of sta¬ 
ble and dense samples of reactive ultracold molecules is 
the ability to understand and control of the loss pro¬ 
cess. It was predicted m. that reactive collisions of 
polar molecules can be stabilized by confining them in 
a quasi-two-dimensional (quasi-2D) geometry with the 
dipoles aligned perpendicular to the 2D plane. Reduc¬ 
tion of reactive collisions is explained by the fact, that 
the molecules are ordered in a way, that they repeal each 
other, hence the probability of a short-range collision is 
significantly reduced. This reduction has been confirmed 
by theoretical calculations based on adiabatic poten¬ 
tials [20) . multichannel quantum dynamics for fermionic 
KRb priES] and by analytical models for universal po¬ 
lar molecules with unit reaction probability at short dis¬ 
tances [22]. A similar effect can be achieved by confin¬ 
ing polar molecules in a quasi-one-dimensional (quasi- 
ID) trap with dipole moments oriented in the direction 


perpendicular to the trap axis. This however does not 
lead to any further reduction of reaction rates in com¬ 
parison to quasi-2D case [24) . 

In this paper we characterize elastic and reactive col¬ 
lisions in quasi-2D traps both for universal and nonuni¬ 
versal polar molecules. The latter have probability of a 
short-range reaction smaller than one, which is described 
by a single dimensionless parameter. Hence, our results 
are applicable to a broad range of polar molecules. We 
investigate how the collision rates change with tuning 
of the characteristic trap length, dipolar length and the 
collision energy. Considering heteronuclear bialkali po¬ 
lar molecules that are at the heart of experimental in¬ 
vestigations, KRb, LiNa, LiK, LiRb and LiCs have al¬ 
lowed reaction channels [23, henceforth their collisions 
should be universal. On the other hand, NaK, NaRb, 
NaCs, KCs, and RbCs do not have reactive channels, so 
they are expected to exhibit nonuniversal collisions in 
their ground state. In general, even nonuniversal polar 
molecules or homonuclear dimers can exhibit universal 
properties when prepared in an excited rotational or vi¬ 
brational state [23 [22]. While universal molecules do 
not have scattering resonances, because their close to 
threshold bound states vanish too fast, nonuniversal sys¬ 
tems should exhibit dense spectrum of overlapping reso¬ 
nances. The latter case can be treated by models comb¬ 
ing quantum-defect methods with assumptions derived 
from random-matrix theory [23 US] • 

The control of reactive collisions has been recently 
demonstrated for fermionic '^°K®^Rb molecule in the rovi¬ 
brational ground state [H El 130]. In the absence of 
an external electric field, molecular interaction is dom¬ 
inated by isotropic van der Waals potential, and ultra¬ 
cold fermionic molecules can either collide is p or s wave 
depending if they are prepared in identical or noniden¬ 
tical internal spin state, respectively. As the ultracold 
p-wave collisions require tunnelling through a centrifugal 
barrier, the identical fermions reactions are significantly 
reduced at ultralow temperatures in comparison to dis- 
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tinguishable particles HU. For a sample polarized by 
some external electric field, the reaction rates in three 
dimensional traps were increasing with the value of the 
induced dipole moment m- Confining molecules in ID 
deep optical lattice creates an effective array of quasi-2D 
traps [30] • In this case polarizing the molecules in the 
direction parallel to the axis of quasi-2D traps, reduces 
the reactive collisions in accordance with the theoreti¬ 
cal models, provided the molecules were prepared in the 
identical spin states [30] ■ On the contrary, for collisions 
in two distinct spin states when the s-wave collision chan¬ 
nel is available, inclusion of an external electric field only 
increases the reactive rates similarly to three-dimensional 
case. 

This work is structured as follows. In section II we 
present our theoretical model that is used to describe re¬ 
active collisions in quasi-2D traps. In particular, we dis¬ 
cuss the limit of a zero-dipole moment when the collision 
rates can be determined analytically from the properties 
of the van der Waals potential. In section III we present 
numerical results based on mutichannel wave function 
propagation. We discuss dependence of the collision rates 
on the dipolar length, confinement length and the kinetic 
energy. Finally we conclude in Sec. IV. 


II. THEORETICAL BACKGROUND 


Interaction of polar molecules differs significantly from 
the interaction of neutral atoms. In case of neutral atoms 
short-range interaction dominates and isotropic van der 
Waals potential is characterized by a Cq coefficient. On 
the other hand polar molecules interact via a long-range 
and anisotropic potential (it depends on the orientation 
of molecules in space) 

(I) 

where d{F) is a dipole moment induced by an external 
electric field F and 9 is an angle between dipoles. In order 
to calculate the collisonal dynamics of polar molecules we 
have used quantum defect method [13 [18] . It assumes 
the separation of lengthscales in the described problem, 
that is visualized in Fig. [^ 

Elastic collision rates {JCf^ and respectively) can 
be described by a complex scattering length El: which 
is directly connected to the scattering matrix S 

di^iE) = aem{E)-iPtm{E) = (2) 


where S£rn,em are diagonal elements of the scattering ma¬ 
trix, im are angular-momentum quantum numbers, E is 



FIG. 1. Lengthscales separation in two polar molecules colli¬ 
sion. Ro, Rioas - the lengthscale of chemical reactions (losses), 
Rn - van der Waals length, ah - characteristic length of the 
trap, ak - de Broglie wavelength. 


a collision energy. Henceforth collision rates are given by 

|1 - 5,™.£m(i?)|' = 2g—|h,„(fc)|V^m(fc) , 

fi 

(3) 

ETm{E) =9^{l- \SimME)?) = - 

IJj 

(4) 


where /r is a reduced mass, k denote the wave vector, g = 
1 besides the case, when two particles are identical {g = 
2). Angular momentum is restricted to even numbers in 
case of identical bosons and to odd numbers in case of 
fermions in an identical state. Moreover fimik) is given 
by a formula 




1 

1 -I- k'^\dem{k)\'^ + 2kPem{k) 


(5) 


For low energies fim{k) —>■ I, when fc —)■ 0, given 
that k\dimik)\ <C 1 and kjStmik) ^ 1. In general 
0 < fem(k) < 1. On the ground of above equations 
(§ and Q one can state that at ultralow temperatures 
imaginary part of the scattering length is responsible for 
reactive scattering rates and modulus of the scattering 
length is responsible for elastic scattering rates. 

We note that for small molecule separations wave func¬ 
tions can be written in a WKB like form 


^'(r 


exp[—i k{x)dx — tip] 

'/(Hr)) 

1 — y exp[f H k{x)dx -I- i(p] 


( 6 ) 


l + y 


/(Hr)) 


where y gives an amplitude of the reflected wave func¬ 
tion. For y = 1 there is no reflected probability flux and 
short range reaction has a unit probability. Whereas for 
?/ = 0 all incoming particles are reflected. At short dis¬ 
tances van der Waals interaction is dominating, so we 
can use well known solutions for the zero energy radial 
wave function in a power-law potential |31j 

^ -|-tan((/?)V2f+i/4 

( 7 ) 



/rRi^E=o cx J2e+i/4: ^ 
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TABLE I. Coefficients pre sent in the formula for universal 
reactive scattering rate ( |10[ ) for different dimensionalities and 
statistics. 


where Rq = (2iJ,Cejh?Y^^. Hence dimensionless scatter¬ 
ing length s = a/d is connected to the short range phase 
by 


/T^ cosjip - 7r/8) 
sm{(f + 7r/8) ’ 


( 8 ) 


where we have used the definition of the average scatter¬ 
ing length [32] 


a = - kRs- (9) 

r(i)^ 

It is worth noting, that as was derived in [21], collision 
rates for y = 1 and without the application of the elec¬ 
tric field have similar magnitude in different number of 
dimensions, provided equivalent gas density (cf. tab. |T|. 

( 10 ) 

where j = 0 {j = 1) stands for even (odd) partial waves, 
g depends on dimensionality, is an equivalent three- 

dimensional gas density taking into account at (squeez¬ 
ing in 3 —iV space dimensions), and Lj{K) is a parameter 
that has dimension of a length, and is defined in Table |T] 


III. SCATTERING LENGTH RESONANCES 
AND COLLISION RATES 

In this section we would like to present results of our 
numerical calculations. The multichannel wave function 
was propagated using the Numerov method, where at 
short range we imposed semiclassical boundary condi¬ 
tions Eq. (|^, defined in terms of s and y parameters. We 
focus here on fermions, mostly for comparison with ex¬ 
perimental results for KRb [30] . In order to generalize the 
outcomes for a whole class of elements and molecules all 
quantities are expressed in dimensionless units. Besides 
the units defined before, we have also used the dipolar 
length ttd = fJLcP/{h?) (fJ- - B, reduced mass, d - an in¬ 
duced dipole moment) and a characterstic length associ¬ 
ated with the trap Oh = (H - trap frequency). 

Moreover in the subsequent equations q will denote a 
wave vector of colliding particles. 

Let us start from the analysis of the complex scattering 
length ([^ . For dipolar length of the order of the trapping 



FIG. 2. The real part of the scattering length for different 
short-range reaction probabilities and the nonzero induced 
dipole moment for qa — 0.05. 


length, the amplitude of its real part is much larger than 
the amplitude of its imaginary part (it can be observed on 
Figs. and [^. Hence, with help of Eq. ^ we can infer, 
that the real part of the scattering length for ultralow 
energies (here of the order of qa = 0.05, that would be 
equivalent to 50nK for KRb) is responsible for the elastic 
scattering rate, whereas the imaginary part governs the 
reactive collisions. 

Fig. [^depicts the real part of the scattering length for 
the fixed dipolar length and the trap length. We have in¬ 
vestigated confinement-induced resonances. For the unit 
probability of reaction at small separation (j/ = 1) reso¬ 
nances are not present, but with diminishing y parameter 
(responsible for the reactivity and the outgoing probabil¬ 
ity flux) the amplitude of resonances increases. 

A similar behavior can be observed for the imaginary 
part, which is shown in Fig. [^ It may be counterintu¬ 
itive that with lowering the reaction probability at small 
distances, the reactive collision rate may increase close 
to the resonance. It can be easily explained. For small 
y parameter the molecules can stay for a longer time at 
short distances and in case of their energy being in reso¬ 
nance (possible only for y < 1) the collision rate can be 
highly enhanced. What is interesting, a resonance, that 
for pure van der Waals interaction occurs at s = 2 for 
fermions, has been strongly shifted by the tight trapping 
potential and the applied electric field to s « 0.4. 

Let us compare the previous outcomes with Fig. [^ 
Here, the imaginary part of the scattering length is de¬ 
picted as a function of s parameter, but also for different 
traps and induced dipole moments. Firstly, we can ob¬ 
serve, that the resonance is shifting not only with reac¬ 
tivity, as was seen before, but also with growing induced 
dipole moment towards lower s values. Moreover, for the 
strong enough applied electric field, the resonance is al¬ 
most totally damped. Similarly with growing squeezing 
in z direction, the amplitude of the resonance gets lower 
and shifted (mind the scale on y axis in Fig. |^. Hence, 
we were able to qualitatively verify, that for sufficiently 
squeezed induced dipoles the head to tail [30] collisions 
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FIG. 3. The imaginary part of the scattering length for differ¬ 
ent short-range reaction probabilities and the nonzero induced 
dipole moment for qa = 0.05. 


are being inhbited, because they are energetically not 
favourable and side by side collisions do dominate. The 
limitation of resonances may be partly understood by the 
fact that in side by side collisions the interaction potential 
is highly repulsive and molecules do not have a chance to 
meet at short distances and react. The outcomes without 
the electric field are compared with an analytic approx¬ 
imation (short dashes on Fig. for a Dirac delta type 
interaction [33) . 

In Fig. we show parametric plots of the real and 
imaginary parts of the scattering length as a function of 
the phase parameter s, for a fixed value of y. It turns 
out that for not too strong values of the dipole moments, 
they are circles with radius depending on reactivity at 
short distances. For universal systems {y = I) the cir¬ 
cle reduces to a single point, whereas with a diminishing 
reactivity the radius increases. It means that for differ¬ 
ent short-range phases (so also s parameters) for small 
short-range reactivity large reactive and elastic rates may 
appear. Bigger dipolar length leads to shifting of the cir¬ 
cles towards larger values of the real part of the scattering 
length and also smaller values of the imaginary part of 
the scattering length. Meanwhile the circle shrinks. We 
can conclude that with help of equation for reactive and 
elastic collision rates ^ even without the precise knowl¬ 
edge of s parameter, that picks a certain point from the 
circle for a certain short-range reactivity it is easy to find 
a range of possible reactive and elastic collision rates. 
What is more for large applied electric helds the system 
approaches the universal limit immune to the short-range 
parameters. 

On Fig. we have shown collision rates for different 
values of s and y parameters as a function of dipolar 
length. Elastic collisions despite differences in the scat¬ 
tering length, approach the universal limit (y = 1) for 
sufficiently high external electric fields, so they become s 
parameter independent. Similarly reactive collisions for 
sufficiently high dipolar length do oscillate around the 
universal curve. The oscillations that are present are the 




FIG. 4. Imaginary part of the complex scattering length for 
different induced dipole moments and the trapping strength. 
Upper panel: an/a = 1.7, lower panel: ah/a = 5.2, short 
dashes: analytical approximation (equation (42) in [^) for 
ad/a = 0, continuous line: numerical calculations for adja = 
0, dots: adja = 0.32, dashes: adja = 0.73, dots and dashes: 
adja — 2.02. All calculations for y = 0.7. 


effect of shape resonances. So for strong dipoles elastic 
collisions dominate and it should be possible to stabi¬ 
lize the system and for example conduct the evaporative 
cooling. 

Now, we will look closer at the dimensionality of the 
system. We investigate how strong should be the dipole 
moment to observe a crossover from a 3D to the quasi- 
2D regime. Let us focus on Fig. For a strong trapping 
(blue lines) and weak dipoles reactive collisions are dom¬ 
inating, but for larger dipolar lengths their amplitude is 
diminishing and elastic collisions win. This is completely 
another way for a weak trap with a large characteristic 
lengthscale (red lines). Here, reactive collisions are for 
small applied electric fields less favourable than elastic 
ones, but finally they aquire similar level. This growth is 
connected to the fact that there is a space for molecules 
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FIG. 5. Parametric plots of the imaginary and the real part 
of the scattering length for all possible values of s parameter. 
Note a different scale on x and y axes in the bottom right 
panel. 



FIG. 6. Reactive and elastic collision rates as a function of the 
ratio of dipolar length to the characteristic trap length and 
the phase given by the s parameter for intermediate values of 
experimental parameters. 


to collide head to tail, that enhances reactivity. The the¬ 
oretical outcomes presented in this paper are valid for a 
broad class of dipolar molecules. Similar conclusions may 
be inferred from experimental investigations m, where 
the presence of the tight trapping potential in the plane 
perpendicular to the external electric field significantly 
supresses the reactive to elastic rate ratio. 

Finally we will check how reactive collision behave as 
a function of the kinetic energy of the colliding particles. 
It is depicted in an dimensionless units in Fig.[^ For two 
different dipolar lengths the dependence is linear, in an 



FIG. 7. Reactive and elastic collision rates as a function of the 
ratio of dipolar length to the characteristic trap length and 
the geometry of the trap for trap parameters corresponding 
to a very weak trap (a/uh = 0.02, dash-dotted and dotted red 
lines for reactive and elastic collision rates, respectively) and 
for moderate trap parameters (a/ah = 0.2, solid and dashed 
blue lines for reactive and elastic collision rates, respectively). 
The figure presents the influence of the trap dimensionality 
on the elastic and reactive rates. 



(aq)" 


FIG. 8. Reactive collision rates as a function of the scatter¬ 
ing kinetic energy. Linear dependence is clearly visible for 
nonzero dipole moments. 


agreement with a Wigner law for an inelastic collisions 
for £ = 1 and with theoretical predictions for universal 
collisions IM] . 


IV. CONCLUDING REMARKS 

We have investigated reactive and elastic collisions in a 
quasi-2D geometry in the presence of an external electric 
field. We have studied the influence of the trapping and 
the induced dipole moment on the resonances of the com¬ 
plex scattering length and collision rates. We have shown 
that for sufficiently large induced dipole moments elastic 
collisions dominate over reactive ones and the system is 
stabilized. Furthermore, we have shown, that reactive 
collision rates scale linearly as a function of the scat- 
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tering energy even for nonzero induced dipole moments. 
For moderate and small dipole moments parametric plots 
of real and imaginary parts of the scattering length for 
various short-range phases are circles with diameter di¬ 
minishing with increasing dipole moment. This is similar 
to calculations without the external electric field [M] and 
shows that even without the knowledge of the short-range 
physics we could predict the possible ranges of elastic and 
reactive scattering rates and their ratio. 
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